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Abstract. 

The Cepheid period-luminosity (P-L) relation is regarded as a linear rela- 
tion (in log[P]) for a wide period range from ~ 2 to ~ 100 days. However, several 
recent controversial works have suggested that the P-L relation derived from the 
Large Magellanic Cloud (LMC) Cepheids exhibits a non-linear feature with a 
break period around 10 days. Here we review the evidence for linear/non-linear 
P-L relations from optical to near infrared bands. We offer a possible theoreti- 
cal explanation to account for the nonlinear P-L relation from the idea of stellar 
photosphere - hydrogen ionization front interaction. 



1. Introduction 



The Cepheid period-luminosity (P-L) relation takes the form of M\ = a\ log(P)+ 
b\, where a and b are the slope and zero-point (ZP) for the P-L relation at a given 
A bandpass, respectively, and P is the pulsation period in days. The Cepheid P- 
L relation is the first rung of the distance ladder and has been widely applied to 
the distance scale studies, which d eliver less than 10% accuracy in distance (for 
example, see lFreedman et al. 2001 ). The Cepheid P-L relation can also be used 
in the stellar pulsation and evolution studies, by comparing the theoretical P-L 
relations to empirical results. The Cepheid P-L relation has been assumed to be 
linear and universal for a long time. Hence Cepheids in the Large Magellanic 
Cloud (LMC) were used to derive the P-L relation due to its nearby proximity 
and the prese nce of a large number equal-distant Cepheids (for example, see 
Udalski et aDll999l b 

Nearly a century after its discovery by Leavitt in 1912, there are still 
some open questions regarding the Cepheid P-L relation. One of them be- 
ing the metallicity dependency of the ZP. The latest empirical d etermination 
of th is dependency is: 7 = -0.29 ± 0.09(12) ± 0.05(S)mag dex" 1 (jMacri et al. 
2006). Furthermore, recent works have s trongly suggested that the Cepheid 



P-L relation could be neither universal (iTammann. Sandage Reindl 2003: 
Ngeow k Kanbur|[2003:lFiorentino et al.ll2007f) nor linear dKanbur Ngeowll2004 
20061: [Sandage. Tammann fc Reindll 12004 : iNgeow et al.1 120051 : iNgeow fc Kanburl 
2006al ). in contradiction to the usual assumptions for the Cepheid P-L rela- 
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tion. In this Proceeding, we concentrate on the discussion of the non-linear P-L 
relations as observed from the LMC Cepheids. 



2. The Non-Linear LMC Cepheid P-L Relation 



Even though almost all of the distance scale studies assume that the LMC P-L 
relation is linear, high quality data for a large number (~ 10 3 ) of fundamental- 
mode LMC Cepheids have provided evidence to imply that the optical LMC 
P-L relations are not linear, with a discontinuity in the slope seen around P = 
10 days (see references given in the Introduction). The justification and the 
significance of selecting the break period at 10 days has been discussed and 
given in the same references above. Various tests have been performed to look for 
the causes of this non-linearity, which include observing strategies, photometric 
errors, extinction errors, influence of outliers, number of long period Cepheids in 
the samples and contamination by overtone Cepheids. However, none of these 
or any combination of the m is found to be responsible f or the observed non- 
linear LMC P-L r elation (iKanbur Ngeowl |2004 120061 : ISandage et all 120041 : 
Ngeow et al.1 120051 : iNgeow fc Kanburl l2006bf ) . This is in contradiction to the 



existing paradigm that the Cepheid P-L relation should be linear. 

Due to intrinsic dispersion along the P-L relation (roughly ~ 0.2mag in 
V band) caused by the finite width of the instability strip, the non-linearity 
of the P-L relatio n is difficult to vi s ualize. This is convincingly demonstrated 
from Figure 1 of[Ng eow Kanburl (|2006bl ). which presents two simulated P- 
L relations constructed from an intrinsicly linear and non- linear P-L relations, 
respectively. Therefore, careful statistical an alyses are required to de tect the 
presence of non-linearity in the P-L relation (jNgeow Kanburl [2006bl ^ . These 
have included the -F-test, non-parametric methods such as LOESS, regression 
methods robust to outliers, the testimator approach and Bayesian Information 
Criterion methods. All have supported the existence of the non-linear LMC 
P-L relat ion at more than a 95% confidence level and a "break" period around 
10 davs (iKanbur fc Ngeowl 12004 120061 : INgeow et al J 120051 : iKanbur etHI 120071 : 
Ngeow et al.ll2008l ). 



I n Table 1, we summarize the F -test results (taken from IKanbur &z Ngeow 
2004 120061 : INgeow et al1l2005l . 120081 : INgeow & Kanburll2008r ) in different band- 



passes from various datasets. The fact that the non-linear results were obtained 
from two different datasets, namely the OGLE (Optical Gravitational Lensing 
Experiment) and the MACHO (MAssive Compact Halo Objects) data, rein- 
force the existence of non-linear LMC P-L relations. Simple black-body argu- 
ments with Ce pheid-like temperatures c an be used to explain the results given 
m Table 1 (see H geow Kanburl [2006bl . for more details). This suggests that 
Cepheid temperature (or equivalently the color) plays an important role in the 
observed non-linear P-L relation. In addition to the P-L relation, the period- 
color (P-C) relation and the instab ility strip were also found to be non-linear in 
the optical for the LMC Cepheids (|Kanbur Ngeowl [2003 ; ISandage et al.ll2004 
Kanbur fc Ngeowl l2006fi . 
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Table 1. Summary of the F-test results. 



Dataset Bandpass Result 



OGLE 




B 


Non-linear 


OGLE 




V 


Non-linear 


MACHO 




V 


Non-linear 


MACHO 




R 


Non-linear 


OGLE 




I 


Non-linear 


OGLE + 


2MASS 


J 


Non-linear 


MACHO 


+ 2MASS 


J 


Non-linear 


OGLE + 


2MASS 


H 


Non-linear 


MACHO 


+ 2MASS 


H 


Non-linear 


OGLE + 


2MASS 


K 


(Marginally) Linear 


MACHO 


+ 2MASS 


K 


(Marginally) Linear 


OGLE + 


IRAS 


3.6/im 


Linear 


OGLE + 


IRAS 


4.5/xm 


Linear 


OGLE + 


IRAS 


5.8/im 


Linear 



3. The HIF - Photosphere Interaction 



The hydrogen ionization front (HIF) is a region of rapid temperature change 
(together with a sharp rise of the opacity) near the surface of Cepheids, where 
hydrogen is partially ionized with a characteristic temperature. Due to its radial 
pulsation, the HIF will move in-and-out in mass distribution for a given Cepheid, 
and both of the HIF and photosphere are not co-moving. Therefore, at certain 
phases of the pulsation, the photosphere of the Cepheid (defined as a layer with 
optical depth of 2/3) can be located at the base of HIF. The "opacity wall" from 
the HIF prevents the photosphere moving in further into the mass distribution. 

Depending on the period, metallicity and the density, the photosphere can 
interact with HIF when it is located at the base of HIF, while the density is 
reasonably low. In this situation, the temperature of the photosphere is the 
characteristic temperature of the HIF which ionizes hydrogen. A consequence 
of this HIF-photosphere interaction is that the P-C relation will be flatter (i.e., 
more independent of the pulsation period) for phases at which this interac- 
tion occurs at low densities. The flat P-C relation has been observed for Ion. 



period Galactic and LMC Cepheids at the maximum light (1 Simon et al.lll99 



Kanbur fc N geowll2004l: [Kanbur et alj|2004t IKanbur Ngeow!l2006l) and fo r RR 



Lyrae at the minimum ligh t (jSturch 1966 : Kanburl 1995: IKanbur fePhillipsl [l996 ; 
Kanbur &: Fernandol 12003 ) . Since the P-L and P-C relations are 2-dimensional 
projection of the P-L-C relation, and the P-L relations normally stated in the 
literature (such as those tested in Table 1) are mean light relations (i.e., aver- 
age over the pulsation phases), the "abnormality" of the P-C relation at certain 
phases could lead to the observed non-linear mean light P-L relation. 

To test the effect of HIF - photosphere interaction on the P-C relation, we 
constructed sequenc es of stellar puls ation models appropriat e for the Galactic 
(|Kanbur et al.1 120041 ) and the LMC (IKanbur fc Ngeowl l2006n Cepheids. From 
the models, we identified the locations of the photosphere and HIF in the mass 
distribution, q = log[l — M(r)/M], where M(r) is the mass within radius r and 
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M is the total mass of a given Cepheid model. A qualitative way to indicate 
the HIF-photosphere interaction is by plotting the "distances" between HIF and 
photo sphere, in terms o f q, ag ainst the pulsational periods as shown in Figure 
15 of iKanbur fc Ngeowl (|2006l ). From that Fi gure, the HIF is located near or 
at the base of the HIF at the maximum light for both of the Galactic and 
LMC models, which lead to the observed flat P-C relation at maximum light. 
However at minimum light, photospheres from the Galactic models are clearly 
far away from the HIF in mass distribution, and hence no interaction occur 
with the HIF. The LMC models, on the other hand, show that the photosphere 
is disengaged from the HIF at minimum light for models with period greater 
than 10 days. In contrast, the short period LMC models suggest the HIF - 
photosphere interaction occurs at most phases. This different behavior of the 
long and short period LMC models may be responsible for the observed non- 
linear LMC P-L relation at mean light. 



4. Conclusion 



Using the F-test, LMC Cepheid P-L relations are found to be non-linear in 
the optical and near infrared JH bands, and become linear around K band 
and longer wavelength, as summarized in Table 1. These results, from a simple 
black-body argument, suggest temperature plays an important role in explaining 
the non-linear P-L relation. One mechanism that can affect the photospheric 
temperature at certain phases of the pulsation is the HIF - photosphere inter- 
action. This interaction can produce a flatter P-C relation at maximum light as 
observed in the Galactic and LMC Cepheids. Our Cepheid models also imply 
that a distinct behavior between the Galactic and LMC models at various phases 
can lead to the observed non-linear LMC P-L relation at mean light. 

The discovery of the non-linearity in the LMC P-L relation has challenged 
the long standing assumption that the Cepheid P-L relation is linear. The 
entire Cepheid based distance scale is based on this assumption and the sit- 
uation needs to be addressed. A more accurate (< 5%) measurement of Hq 
than exists at present is crucial because it will h elp to break the degeneracy be- 
tween Qm and Hq in CMB measurements (e.g.. iTegmark et al"]|2004l ; lHull2005l ; 



Spergel et al.l 120071 ) , required in the era of precision cosmology. The optimal 



way to d o this is to calibra t e the C epheid distance scale to an accuracy of a few 
percent. Ngeow fc Kanbun (2006c) show that the error in estimating Hq arising 
from using a linear P-L relation when the true relation is intrinsically non-linear 
is about 2%. As other sources of larger systematic error in Hq are addressed (e.g. 
the absolute calibration or the metallicity dependence) and a number of recent 
attempts at redu cing zero point errors on the Cepheid distance scale (e.g., see 
Macri et alJ l2006). our lack of knowledge about the slope of the P-L relation will 



become increasingly important. Equally as important, a proper understanding 
of the non- linear P-L relation is vital for stellar evolution/pulsation studies of 
Cepheids, to develop a theory explaining the physics behind these observations. 
One such theory is the proposed HIF - photosp here interaction as described in 
previous section and in IKanbur &; Ngeowl ( 20061 ). 
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